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An Upper Bound on the Flux Ratio of p CrB's Companion at 

1.6 /im 

C. Bender 1 , M. Simon 1 , L. Prato 2 ' 3 , T. Mazeh 4 , and S. Zucker 5 ' 6 

ABSTRACT 

We use high resolution infrared spectroscopy to investigate the 2001 report 
by Gatewood and colleagues that p CrB's candidate extrasolar planet companion 
is really a low-mass star with mass 0.14 ± 0.05 M Q . We do not detect evidence 
of such a companion; the upper bounds on the (companion/primary) flux ratio 
at 1.6 fim are less than 0.0024 and 0.005 at the 90 and 99% confidence levels, re- 
spectively. Using the if -band mass-luminosity relationship calculated by Baraffe 
and colleagues, the corresponding upper limits on the companion mass are 0.11 
and 0.15 M . Our results indicate that the infrared spectroscopic technique can 
detect companions in binaries with flux ratios as low as 0.01 to 0.02. 

Subject headings: binaries: spectroscopic — stars: individual (p Coronae Bore- 
alis) — techniques: radial velocities, spectroscopic 

1. Introduction 

p Coronae Borealis (p CrB, HIP 78459), a well-studied GOV star with an estimated 
mass of 1.0 M and an age of 10 Gyr, was among the first stars identified with an extrasolar 
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planet (ESP) candidate (Noyes et al. 1997, hereafter N97). If the mass of the companion, m, 
is small relative to p CrB, the system's mass function indicates that m sin i = 1.1 Mj upiter , 
where i is the orbital inclination. For a random distribution of orbital plane inclinations, 
the expectation value of m is (4/7r) x 1.1 M Jupiter , justifying the interpretation that p CrB's 
companion is an ESP. 

Subsequently, Gatewood, Han, & Black (2001, hereafter GHB) used data from Hipparcos 
and the ground-based Multichannel Astrometric Photometer (Gatewood 1987) to measure 
the astrometric signature of the reflex motion of p CrB. They derived an astrometric orbit 
with a semi-major axis of 1.65 ± 0.33 mas and an inclination of i ~ 0.5°, indicating that the 
companion is a late type M-dwarf star with mass 0.14 ± 0.05 M Q . In a large sample of ESP 
candidates, only a few systems seen at low inclinations are expected. It was therefore sur- 
prising that, in a companion paper, Han, Black, & Gatewood (2001, hereafter HBG) found, 
by combining Hipparcos data with the spectroscopic orbital elements of 30 ESP candidates, 
that a significant fraction may have small inclinations. The semi-major axes of the orbits 
that HBG derived are, however, at the limits of the precision of Hipparcos measurements. 
Pourbaix & Arenou (2001) and Zucker & Mazeh (2001) demonstrated that in these circum- 
stances the small i's were an artifact of HBG's analysis. Indeed, McGrath et al.'s (2002) 
HST astrometric observation of p 1 Cnc (55 Cnc), a star estimated by HBG to have a partic- 
ularly low inclination orbit, placed an upper limit on its reflex motion, confirming that its 
companion is sub-stellar. 

We have been using high-resolution infrared (IR) spectroscopy to detect low-mass com- 
panions in spectroscopic binaries first identified in visible light (Mazeh et al. 2002, hereafter 
M02; Prato et al. 2002, hereafter P02). If p CrB's companion had the mass of a late M 
spectral type main-sequence star, we expected it to lie within the sensitivity limits of the 
technique (e.g. Mazeh et al. 2003). We describe here our search for the companion in the 
p CrB system by IR spectroscopy. §2 details extension of the IR technique to p CrB, our 
observations, and data reduction. In §3 we describe our analysis and the measurements of 
p CrB. In §4 we set an upper bound on the 1.6 pm flux of p CrB's companion. We calcu- 
late an upper limit on the companion mass using theoretical models and comment on the 
significance of our result in §5. 



2. Observations and Data Reduction 

We observed p CrB during four runs in the spring of 2001 (see Table 1) with the W. 
M. Keck Observatory's cross-dispersed echelle IR spectrometer NIRSPEC (McLean et al. 
1998, 2000) on the Keck II telescope with and without adaptive optics (AO). We used 
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NIRSPEC in its high resolution mode centered at 1.555 /im. This provided 9 spectral orders, 
45 - 53, from 1.450 - 1.705 //m with resolutions of 28,000 and 35,000 in non-AO and AO 
modes respectively, for a 2 pixel slit width. Further details of the instrument setup and 
observation procedure are given in P02. We extracted the spectra using the RED SPEC 1 
software package, and determined the dispersion solution using OH emission lines of the 
night sky, recorded simultaneously with the observations and identified from the catalog of 
Rousselot et al. (2000), and arc lamp spectra. 

Application of the IR spectroscopic technique to the detection of low mass secondaries 
requires good templates for the primary and the secondary. Observations of these templates 
took place during 2000 and 2001 (P02). Prior to this work (e.g., P02, M02), we used spectra 
in only NIRSPEC order 49 because it is nearly completely free of terrestrial absorption. 
However, the sensitivity of our analysis depends on the number of stellar spectral lines 
available and hence on the spectral range of the data. To reach the low flux levels expected 
for p CrB's companion we used additional spectral orders. Figure 1 shows the spectrum of 
the AO star HD 13372 in orders 45 - 53. The stellar spectrum contains only weak, broad 
lines of the hydrogen Brackett series; all other absorption shown in the figure is terrestrial. 
The absence of terrestrial absorption in order 49 is apparent. Order 46 contains terrestrial 
absorption in the P and Q branches of a CH 4 band; orders 47 and 48 contain C0 2 absorption 
lines. Orders 51 - 53, and to a lesser extent order 50, are contaminated by many terrestrial 
absorption lines. In order 45, with our choice of central wavelength, starlight illuminated the 
array in only one of the two telescope nod positions. We concentrated, therefore, on trying 
to correct atmospheric absorption only in the spectra of orders 46, 47, and 48. 

Terrestrial absorption lines are usually removed by dividing the target spectrum by a 
featureless stellar spectrum observed nearly simultaneously at the same airmass. We took 
a different approach because we had to correct many template spectra obtained in different 
observing sessions. We expect that the C0 2 and CH 4 absorption bands are independent of 
weather and vary only with airmass. We therefore obtained NIRSPEC spectra of five AO 
stars on July 17, 2002, over the range of airmasses of our p CrB and template spectra. We 
removed the Brackett lines by fitting them with Lorentzian profiles and used the terrestrial 
absorption lines in orders 46, 47, and 48 to register each observation on a sub-pixel scale. 
For each of the terrestrial lines, we fitted the dependence of the absorption on airmass, and 
were thus able to calculate terrestrial absorption spectra at the exact airmasses of the p CrB 
and template spectra. This procedure worked well for order 47 and less well in orders 46 and 
48, probably because their telluric lines are deeper and have more complex structure than 
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those in order 47. The residuals in orders 46 and 48 appear in only a few percent or less 
of free spectral range of these orders so we considered them negligible. Figure 2 shows the 
order 46 - 49 p CrB spectra measured on June 2, 2001, with orders 46, 47, and 48 corrected 
for telluric absorption. 

To ensure consistent results, we reextracted orders 46 - 49 for the entire template set, 
calculated the dispersion solution using the OH lines of the night sky, corrected the spectra 
for terrestrial absorption, and rederived the stellar radial velocities following the technique 
described in P02. Figures 3-6 show the template spectra for orders 46 - 49 in the laboratory 
reference frame, with the terrestrial absorption lines removed from orders 46, 47, and 48. 

3. The Companion/Primary Flux Ratio of p CrB 

Our procedure to search for evidence of p CrB's companion is similar to that described by 
M02 for spectroscopic binaries using the two-dimensional cross-correlation routine TODCOR 
(Zucker & Mazeh 1994). The analysis identifies the templates that provide the best matching 
primary and secondary, the ratio of the secondary and primary fluxes, a, and the radial 
velocities of the components. The measured amplitude of p CrB's reflex velocity is 67 ms" 1 
(N97). If its companion had the 0.14 M mass reported by GHB, the maximum velocity 
difference of the primary and secondary would be 0.4 km s -1 . This is less than the lkms -1 
velocity precision of our technique (P02), so we did not expect to detect a velocity difference. 
To obtain a velocity difference large enough for us to detect reliably, > 1 km s _1 , would require 
that the companion have mass < 0.05 M ; the flux of such a substellar object blended with 
a G star would be too small for us to detect spectroscopically. We therefore regarded the 
velocity difference of the components as fixed at and focused on measuring the system's 
flux ratio. 

We wrote a routine to model binaries from pairs of templates at prescribed flux ratios, to 
cross-correlate them with our p CrB spectra, and to identify the model most closely matching 
p CrB. With this, we analyzed separately the p CrB spectra in orders 46 - 49 for each of the 
four observing sessions using a wide selection of templates. The highest correlation values 
resulted from using HD 4614 or GL 160 for the primary template, depending on the particular 
p CrB spectrum, and the late M-type templates for the secondary. Changing the secondary 
spectral type by a few subclasses resulted in negligible changes in the correlation, and we 
concluded that we could only discriminate between early and late type M-star secondaries. 
Based on this, we averaged the a's resulting from the late type secondaries GL 406, LHS 292, 
GL 644C, LHS 2351, and LHS 2065 for each p CrB spectrum. These averaged values are 
our measured flux ratios and are listed in Table 2. The arithmetic mean of the measured 
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flux ratios in Table 2 is —0.0034. This small value, and the large scatter around it, from 
—0.0260 to 0.0164, suggest that the observations have not detected the companion. The 
primary and secondary templates are not likely to match exactly the effective temperature 
and metallicity of the p CrB primary and companion. The fact that the average flux ratio 
is essentially zero indicates that any such mismatch is not causing a systematic error. 

To interpret the apparent non-detection we must derive an upper bound, at a specified 
level of significance, on the flux ratio. The conventional approach would be to calculate the 
mean of the measured flux ratios in Table 2, to assume they follow a Gaussian distribution 
and derive the standard deviation of the mean, and to use this to place confidence limits on 
an upper bound. The measured flux ratios, however, do not necessarily follow a Gaussian 
distribution, and the number of measurements, 16, is too few to test the assumption reliably. 
We chose, therefore, to take a different approach, using the data itself to model a large number 
of observations, and from these directly measure the probabilities of a detection at a certain 
flux ratio, independent of underlying statistical assumptions. §4 describes this approach. 

4. Estimate of the Flux Ratio Upper Bound 

4.1. Model Binaries 

In this section, we estimate the upper bound on the ratio of the companion's flux to that 
of p CrB that follows from the measurements in Table 2. Our approach starts by analyzing 
sets of model binaries. To ensure the models closely represented the original observations, we 
used the observed p CrB spectra as the basis for the model primaries. We subdivided each 
observed p CrB spectrum at wavelengths where there are no lines, into smaller sections with 
lengths from 10 to 60 A. We reassembled these sections in random sequences to create a set of 
500 unique model primaries for each observation, with absorption lines, noise characteristics, 
and free spectral range identical to the original p CrB spectrum. Figure 7 shows examples 
of reassembled model primaries derived from the January 2001 p CrB spectra. We used the 
M7 spectral type star LHS 2351 for the secondary in all of the models, and combined it, 
unmodified, with the randomized primaries at eight flux ratio values, from 0.001 - 0.020, to 
create the model binary spectra. This provided us with 128 sets (4 observations x 4 orders 
x 8 flux ratios) of 500 model binaries each, with which to evaluate the measured flux ratios 
in Table 2. 

We analyzed the model binaries in the same way as the original p CrB spectra (§3). For 
each set of models, we cut and reassembled a matching set of primary templates from either 
HD 4614 or GL 160, depending on which was used in the corresponding p CrB analysis. We 
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restricted the analysis to only the M6.5 spectral type star LHS 292 for the secondary template 
because the correlation value is insensitive to small changes in the secondary spectral type 
(§3). 

The analysis produced 500 "measured" flux ratios for each of the 128 sets of models. For 
a given set of models, the measured flux ratios scattered around the input value. Figure 8 
shows two of these distributions from January 2001, order 49, a = 0.010, 0.005, along with 
the corresponding measured value of —0.0046 from Table 2. 



4.2. Flux Ratio Upper Bound 

From each of our model distributions we can determine the probability that the corre- 
sponding p CrB observation returns a flux ratio as small as the one measured. Consider, 
for example, the results for January 2001, order 49. If the actual flux ratio were 0.010, the 
probability that the observation and its analysis produces a flux ratio as small as the one 
measured, —0.0046, is given by the number of model binaries that returned flux ratios of 
equal value or smaller. Figure 8a shows the model distribution derived for this case; 24 
out of 500 models (probability=0.048) produced an a smaller than —0.0046. If, then, the 
true a were 0.010, the probability of a measurement producing a value as low as or lower 
than —0.0046 is unlikely, and we could consider 0.010 an upper bound with 95% confidence 
(1 — 0.048 ~ 0.95). Suppose that instead the true a were smaller, 0.005. In this case, a 
larger number of models, now 76 out of 500, or 15%, give an a as small as the measured one 
(Figure 8b). We would therefore consider 0.005 as an upper bound with only 85% confidence. 

For each of the 16 modeled observations, we calculated, at the eight input model flux 
ratios, the probability of measuring a flux ratio as low as or lower than the measured value 
in Table 2. To combine those probability values into one value, we use Fisher's method to 
combine independent p- values (Fisher 1932). In this method, one first calculates the product 
of probabilities. This value, according to Fisher (1932), is distributed in a known way, and 
we can use its known distribution to infer the new probability. The product value, k, is 

n-l 

k = Y[ Pl (i) 

i=0 

Fisher (1932) showed that the statistic F = —2 In A; has a x 2 distribution with In degrees 
of freedom. Following Fisher's recipe 2 , we use this known distribution to calculate P, the 
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combined p-value from the 16 probabilities. Figure 9 shows P as a function of input flux 
ratio. The combined probability curve intersects the 99% confidence level at a = 0.005. 
Similarly, the flux ratio 0.0024 can be ruled out with 90% confidence. 

5. Discussion 

We consider the flux ratio from 1.55 to 1.65 /im as equivalent to that in the if -band 
because spectra of late M dwarfs between 1.4 and 1.7 pm are not marked by the very deep 
H2O absorption characteristic of L spectral type dwarfs (e.g. Leggett et al. 2001). The 
2MASS ii-band magnitude of p CrB is 3.99. This is in excellent agreement with its Hipparcos 
V magnitude, 5.39 (HIP 78459), and the V — H color, 1.36, for a GOV spectral type star 
(Tokunaga 2000). With the Hipparcos distance of 17.3 ± 0.2 pc, p CrB's absolute H-band 
magnitude is Mh = 2.8. The flux ratio upper bound (§4.2) therefore corresponds to an 
if-band brightness limit for the companion of M H = 8.6 at the 99% confidence level and 
M H = 9.3 at the 90% confidence level. 

Figure 10 shows the absolute H magnitude versus mass, 10 Gyr isochrone calculated 
by Baraffe et al. (1998) for low mass stars and brown dwarfs, with the model values for 
masses less than 0.1 M Q updated from Baraffe et al. (2003). At the scale of the figure, the 
5 Gyr isochrone is indistinguishable from the one for 10 Gyr. We also plot the locations 
of main-sequence dwarfs whose masses have been measured dynamically; the references are 
cited in the figure caption. For reference, the figure shows the (companion/primary) flux 
ratio upper bounds calculated in §4. The flux ratio upper bound at the 99% confidence 
level, a = 0.005, falls at M ~ 0.15 M Q , close to GHB's reported mass of the companion, 
0.14 ± 0.05 M Q . Similarly, a = 0.0024 sets an upper bound on the mass of M ~ 0.11 M 
with 90% confidence. While these mass upper bounds do not provide either a definitive 
confirmation or rejection of GHB's reported mass of the companion, they do suggest strongly 
that its mass is smaller than their value. 

The analysis described in §3 and §4, and illustrated by the histograms in Figure 8 shows 
that, by increasing the spectral range through the use of several orders, the IR spectroscopic 
technique can detect binaries with 1.6 pm flux ratios in the range of 0.01 to 0.02. This is a 
significant advance over our previous detections at flux ratios of 0.04 to 0.05 (Mazeh et al. 
2003). 
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6. Summary 

1) Our measurements set upper limits of 0.0024 and 0.005 on the ratio of the 1.6 ^m 
flux of p CrB's companion to p CrB, at the 90% and 99% confidence levels, respectively. 

2) Using Baraffe et al.'s (1998, 2003) calculations of the iJ-band mass-luminosity rela- 
tion, these flux limits correspond to mass upper bounds of 0.11 and 0.15 M Q at the 90 and 
99% confidence limits, respectively. 

3) Our analysis of the model binaries tested the sensitivity of our technique when using 
multiple orders. The results indicate that binaries with 1.6 /xm flux ratios as low as 0.01 to 
0.02 are detectable by high spectral resolution IR spectroscopy. 
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Fig. 1. — NIRSPEC spectrum of the AO star HD133772 for orders 45 - 53. The spectra 
have been flattened and normalized to their continuum levels. The stellar spectra are nearly 
featureless except for the very broad hydrogen Brackett lines indicated. The remaining 
absorption lines are terrestrial. A CH 4 band is indicated in order 46 and two CO2 bands are 
marked in orders 47 and 48. 
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Fig. 2.— NIRSPEC spectrum of p CrB from June 2, 2001 for orders 46 - 49. The spectra 
have been flattened and normalized to their continuum levels. The terrestrial absorption 
lines in orders 46, 47, and 48 have been removed by ratioing with calculated spectra (see 
text). 
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Fig. 3. — NIRSPEC order 46 templates. The spectra have been flattened, normalized to 
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Fig. 4. — Same as Figure 3, 



but for NIRSPEC order 47. 
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Fig. 5. — Same as Figure 3, 



but for NIRSPEC order 48. 
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Fig. 6. — Same as Figure 3, 



but for NIRSPEC order 49. 
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Fig. 7. — Example model binary primary spectra from the January 2001 model set. The 
p CrB spectra were subdivided, at wavelengths with no spectral lines, into smaller sections 
with lengths 10 - 60 A. The sections were reassembled randomly to create unique spectra 
with absorption lines, noise characteristics, and spectral range identical to the original p CrB 
spectra. 
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Fig. 8. — a: The distribution of derived flux ratios in a sample of 500 model binaries 
having the characteristics of the p CrB order 49 spectrum measured in January 2001 and 
the companion GL 406 at flux ratio a = 0.010. The dashed line indicates the input a and 
the arrow the measured value, —0.0046 (Table 2). b: Same as (a) but for input a = 0.005. 
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Fig. 9. — The curve shows the combined probabilities that our 16 measurements (4 spectral 
orders on 4 occasions) have a flux ratio less than or equal to the measured value, as a function 
of the model binary input flux ratio. The straight lines with constant probability show the 
90% and 99% confidence levels and intersect the curve of combined probabilities at a flux 
ratios of 0.0024 and 0.005. 
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Fig. 10. — Absolute H-band magnitude vs. mass; the (companion/primary) flux ratio at 
1.6 jum, a, is also shown for reference. The dotted curve corresponds to BCAH98's 10 Gyr 
isochrone, with values for M/M Q < 0.1 updated by Baraffe et al. (2003). The vertical lines 
with constant mass represent GHB's result and la error estimate. The horizontal lines at 
Mh = 8.6 and 9.3 represent our flux ratio limits of 0.005 and 0.0024 for p CrB's companion 
at the 99% and 90% confidence limits, and corresponds to a mass upper bounds of 0.15 and 
0.11 M , respectively. The data points are compiled from the literature and the symbols 
correspond to their source as follows: A - Henry et al. (1999); V - Henry & McCarthy 
(1993); □ - Segransan et al. (2000); x - Forveille et al. (1999); > - Torres et al. (1999); < - 
Martin et al. (1998). M H values come from these and Delfosse et al. (2000), Leggett et al. 
(2002), and Leggett (1992). 



Table 1. Summary of p CrB Observations 



UT Date Obs Mode Total Exp. Time (s) 



2001 Jan 7 AO 640 

2001 Feb 2 non-AO 960 

2001 May 3 non-AO 840 

2001 Jun 2 non-AO 840 
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Table 2. Companion/Primary Flux Ratios at 1.6 fim 



Observation 


TV T T T~\ (~i T "A T — \ f" i 

NIRSPEC 




Date 


Order 


a 


2001 Jan 7 


46 


-0.0100 




47 


0.0125 




48 


-0.0090 




49 


-0.0046 


2001 Feb 2 


46 


-0.0084 




47 


0.0068 




48 


-0.0114 




49 


0.0164 


2001 May 3 


46 


-0.0148 




47 


0.0065 




48 


-0.0148 




49 


-0.0012 


2001 Jun 2 


46 


-0.0260 




47 


0.0118 




48 


-0.0020 




49 


-0.0054 



